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ABSTRACT 

We present a useful analytic approximation to the solution of the Lane-Emden 
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equation for infinite polytropic index - the isothermal sphere. The optimized 
expression obtained for the density profile can be accurate to within 0.04% 
within 5 core-radii and to 0.1% within 10 core-radii. 
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•rn . 1 INTRODUCTION 

We construct an analytic approximation to the full non-singular isothermal sphere. The 
approximations currently in use provide a good fit either to the inner regions (r < 2 r corc ) 
or to the asymptotic behavior. The Lane-Emden equation for the gaseous polytrope with 
polytropic index n — > oo is identical to that of a self-gravitating isothermal sphere. All 
polytropes with n > 5 are infinite and hence no analytic solutions exist. In a recent paper, 
[Liu (19961) has exhaustively examined approximate analytic solutions for polytropes with 
general index, where he obtains a solution for the isothermal case to within < 1%. In this 
brief note, we report a simpler analytic form, accurate to within 0.04% within 5 core-radii. 



2 THE ISOTHERMAL SPHERE 



The Lane-Emden equation written in terms of the standard variables (see Emden 1907| ) is 
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where w = In po/p; £ = ?"/ro; r o = ^V^Gpo; a is the constant velocity dispersion and r 
the core radius. 

The solution for the density distribution can be expanded into a series, 

- = i - (k 2 + (ix 4 - - (2) 

Po 6 45 

The approximation that is used often in the literature (which we denote by the subscript 
(a) henceforth) follows from truncating equation (2) at the second order. 

P(0- = T^T- (3) 
+ 6 

While the mathematical form is simple, this profile overestimates mass outside ro and differs 
from the exact solution by a factor of 3 as r — > oo. 

As an ansatz, we attempt the following analytic form for the approximation: 
A B 

PlW approx — a 2 + £2 £2 _(_ £2 ' W 

where A, B, a 2 and 6 2 are to be determined. The solution obtained (see Appendix for details) 
is given by, 

p(C) b = - — ^ (5) 

'10 ~ 12 

This expression is correct asymptotically and is accurate to within 1% up to £ = 5. Pre- 
serving the general form, this approximation can be optimized further yielding, 

a 2 = [lO-(^-^e)] ; 6 2 = [12 - ( ±5 + ^c)] ; 

4 = (5 + c) ; | = (4 + e), (6) 



where the additional parameters e and 5 have now been chosen (see Appendix for details) 
to optimize the degree of accuracy in terms of agreement with the full exact solution. This 
solution is plotted in Figures 1 and 2. The mathematical form of the optimized solution is 
convenient appropos Abel inversion since it can be inverted analytically. The corresponding 
projected quantities for our formula - the surface density and mass are easily computed and 
have the following analytic forms: 

E(R) = Alx - Bn ( 7 ) 
V } Va 2 + R 2 Vb 2 + R 2 ' 

where R is the projected radius. These spherical models can be easily generalized to describe 
elliptical mass distributions as well, akin to those proposed by |Kassiola, Kovner, fc Fort 
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Figure 1. VARIOUS APPROXIMATIONS TO THE ISOTHERMAL SOLUTION: solid curve - exact solution, dotted curve 
- approx(a), dashed curve - approx(b), long-dashed curve - approx(c) 



(1993). 

M 3D (0 = 4n[Aa£ - tan" 1 ^)) - 5 6(| - {tan" 1 ^))], 
a a o b 

The potential on the plane corresponding to the surface density S is, 

(j) 2D = A7r[Va 2 + R 2 -a\nR - aln(a 2 + avV + R 2 )} 
- Bn[Vb 2 + R 2 -blnR- bln(b 2 + bVb 2 + R 2 )}. 



(8) 



(9) 



These projected quantities are of interest in many physical problems - for instance, in the 
context of modelling gravitational lensing observations. The two primary effects producing 
by lensing are the isotropic manginification and the anisotropic shear. The magnification k 
produced by the potential is given by, 
An Bit 



k(R) = K [ 



Va 2 + R 2 Vb 2 + R 2 



(10) 



where k = ^— ; and S crit is the critical surface density given the geometrical configuration 

•^crit 

of the source, lens and the observer. And the induced shear 7 is, 



i(R) = «o[ 
+ 



An 



VR 2 + a 2 R 2 
Btt 2Btt 



VR 2 + b 2 R 2 



{VR 2 + b 2 -b)]. 



(11) 
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Figure 2. DIFFERENTIAL ERRORS: (i) TOP PANEL - density (dotted curve - approx(a), dashed curve - approx(b), long- 
dashed curve - approx(c)) and (ii) LOWER PANEL - the potential (dotted curve - approx (a), dashed curve - approx(b), 
long-dashed curve - approx(c)) (the dotted curve and dashed curve are coincident) 
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Figure 3. (i) LEFT PANEL: MASS ENCLOSED WITHIN r (solid curve - exact solution, dotted curve - approx(a), dashed 
curve - approx(b), long-dashed curve - approx(c)) and (ii) RIGHT PANEL: circular velocity profiles (solid curve - exact solution, 
dotted curve - approx(a), dashed curve - approx(b), long-dashed curve - approx(c)) 

3 CONCLUSIONS 



The analytic approximation presented above is potentially useful in the context of many 
physical problems and is particularly useful since the projected quantities have simple ana- 
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lytic forms. We also point out that within 5 core-radii, while the analytic approximation to 
the isothermal sphere currently in use systematically over-estimates the mass enclosed our 
formula is accurate to within 0.04%. 



4 APPENDIX 

4.1 Fitting to the exact solution 

Starting with our ansatz for the functional form, 
A B 

P(0approx = fl2 + £2 ~~ tf> + £ 2 ' ( 12 ) 

Expanding the above as, 

P(0a PP rox = £2 ) - £2"), (13) 

P(0approx = 4(1 - S + ^ + -) (14) 

(ji (Jj (Jj 

& 21 6 2 + 6 4 + 

Comparing terms to corresponding orders in £ in equations 1 and 2, we obtain the following 
system of equations: 

A B _ A B 1 A B 1 



a 



2 6 2 'a 4 b 4 6 ' a 6 6 6 45 



Requiring the asymptotic behavior as £ — > oo to match up to the full solution, we obtain the 
additional equation to close the above system of 4 simultaneous equations in 4 unknowns. 

A - B = 2, (16) 

we obtain an exact solution!, 

^ = ire - rre- (17) 

x 10 T 12 

4.2 Optimized approximation 

We introduce additional parameters e, 5, x and y and constrain their values to the desired 
degree of accuracy as follows, 

P(0e = 5+ ,l ~ • (18) 

~"~ 10 - xe 12 - ye 
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We require agreement asymptotically as £ — > oo, therefore simplifying the equation above 
for large e, 

(5 + e)(10 - xe) - (4 + e)(12 - ye) = 2, (19) 
ignoring terms of 0(e 2 ), 

* = {\y- \). (20) 

Substituting the above back into equation (16), 
5 + e 4 + e 

1 + ~ 1 + 



p(0c = : , £2 - g2 ■ (21) 

10 -xe 1 + + 



Now expanding the above two terms on the RHS into a series and matching the corresponding 
terms of the same order in £ to 0(e) and matching asymptotically to S degree of accuracy 
we have, 

y = (— + — )• (22) 

For an RMS error of 0.04% to £ = 5, we find e = -0.07 and 5 = -0.269. If we do not treat 
e as small and do not demand the correct asymptotic solution at large £, we can obtain a 
solution with RMS error of 0.1% at £ = 10. This is e = —0.635 and 5 = —0.6 which is used 
in the plots. In this note, we have approximated the full untruncated isothermal sphere and 
have not addressed the problems of truncating it. 
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